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Executive  Summary 

Following  on  from  the  work  on  the  pyrochlore  materials  (A2B2O7),  we  have  extended 
the  contour  map  approach  by  application  to  a  wider  range  of  perovskite  related 
materials  with  the  general  formula  ABO3.  Intrinsic  defect  processes  were  compared 
within  this  extensive  range  of  A3+B3+03  compounds  with  A  cation  radii  varying  from 
Sc3+  to  La3+  and  B  cation  radii  from  Al3+  to  In3+.  A  change  in  preference  between 
cation  antisite,  Schottky  and  oxygen  Frenkel  disorder  are  predicted  as  a  function  of 
cation  radius  and  trends  have  been  explained  in  terms  of  variations  in  lattice  energy 
and  changes  in  crystal  structure.  Of  particular  interest  was  the  unexpected  result  that 
some  defect  processes  mirror  changes  in  crystal  structure  while  others  change 
smoothly  across  crystal  boundaries. 

This  work  fulfils  a  commitment  in  the  original  contract  to  “Identify  compositions  that 
exhibit  complex  defect  behaviour  in  perovskite  oxides.” 
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Introduction 


Due  to  their  many  actual  and  potential  technological  applications,  perovskite  materials 
have  been  investigated  extensively  [1-6].  Initially  much  of  this  prior  work  considered 
2:4  perovskites,  for  example  the  original  perovskite  compositions  Ca2+Ti4+C>3  and  the 
ferroelectric  phase  Ba2+Ti4+C>3  [7,  8].  Nevertheless  ABO3  compositions  in  which  both 
the  A  and  B  cations  assume  3+  formal  charge  states  are  also  of  considerable 
importance  with  applications  as  microwave  dielectrics  [9],  sensors  [10,  11]  and  solid 
oxide  fuel  cells  (SOFC’s)  [1].  In  all  cases  the  application  is  facilitated  by  inducing 
defects  (in  some  cases  via  dopanting)  within  the  perovskite  lattice. 

Here  our  interest  is  focused  on  intrinsically  pure  materials  as  a  first  step  towards  a 
broader  understanding  of  defects  generally.  For  example,  when  the  lattice  is 
disordered,  defects  on  the  cation  sublattice  will  remain  due  to  their  mobiliy  being 
much  lower  than  those  generated  on  the  oxygen  sublattice  (which  are  known  to  be 
mobile  since  these  materials  are  used  as  anodes  in  SOFC’s  due  to  their  high  oxygen 
ion  conductivities).  One  of  the  more  common  defect  processes  on  the  cation  sublattice 
is  the  antisite  reaction,  where  the  A  and  B  cations  are  swapped.  For  materials  such  as 
(Ca,Sr)Ti03,  where  the  charges  on  the  A  and  B  cations  are  not  equivalent,  a 
considerable  charge  redistribution  will  be  introduced  into  the  lattice  with  intense 
associated  strains.  An  additional  factor  contributing  to  the  lattice  strain  is  the 
significant  size  difference  between  the  A  and  B  cations. 

Here  we  extend  the  studies  of  disorder  processes  in  pyrochlores  to  the  perovsites  by 
considering  ABO3  materials  in  which  the  A  and  B  cations  both  assume  formal  3+ 
valence  states.  One  of  the  primary  motivations  for  studying  these  3:3  compounds  is 
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that  cation  disorder  does  not  lead  to  a  charge  imbalance.  It  may  be  that  the 
equivalence  of  the  charge  states  leads  to  lower  disorder  energies;  however,  for  many 
perovskites  we  are  still  faced  with  a  significant  cation  size  difference. 

We  will  consider  an  exceptionally  broad  range  of  ABO3  compositions  with  the  A 
cation  ranging  from  Sc3+  (1.05A)  to  La3+  (1.34A)  and  the  B  cation  ranging  from  Al3+ 
(0.535A)  to  In3+  (0.8A).  Consequently,  at  one  extreme  of  the  compositional  range  the 
A  and  B  cations  are  of  very  similar  size.  Such  A3+B3+C>3  compositions  no  longer 
exhibit  a  perovskite  like  structure,  but  take  on  the  cubic  bixbyite  structure. 

In  this  study,  atomic  scale  computer  simulation  techniques  are  used  to  predict  the 
structures  and  energies  associated  with  defect  formation  in  these  3:3  compounds. 
Those  structures  which  exhibit  larger  energies  for  disorder  processes  will  assume  a 
greater  residual  or  retained  energy  upon  disorder  and  on  this  basis  are  assumed  less 
tolerant  of  such  defects  [12].  Such  an  approach  has  proved  successful  in  earlier  work 
at  predicting  the  relative  structure  boundary  between  pyrochlore  (in  which  the  cations 
are  ordered)  and  fluorite  (in  which  the  cations  are  disordered)  [12],  The  approach  was 
also  successful  in  predicting  the  order  /  disorder  transformation  in  A2HfC>7 
compositions  [13]  and  nonstoichiometry  across  the  entire  A2B2O7  compositional 
range  [14]. 
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Methodology 

The  calculations  presented  here  are  based  upon  a  classical  Bom-like  description  of  an 
ionic  crystal  lattice  [15].  The  interatomic  forces  acting  between  ions  are  resolved  into 
two  terms:  long-range  Coulombic,  which  were  summed  via  the  Ewald  method  [16], 
and  short-range  pair  terms,  which  were  modelled  using  parameterised  Buckingham 
pair  potentials  [17].  The  perfect  lattice  is  defined  by  tessellating  the  unit  cell 
throughout  space  using  periodic  boundary  conditions  as  defined  by  the 
crystallographic  lattice  vectors.  The  lattice  energy  (El)  can  then  be  expressed  as 
follows: 


UL  = 


4  TIE. 


'^+Aevc 


0  i*j 


r>j 


r6. 

V 


Equation  1 


where  Ay,  py  and  Cy  are  the  adjustable  parameters,  ry  are  the  interionic  separation  and 
qi?  qj  are  the  charges  on  ions  i  and  j  respectively.  The  parameters  of  the  short-range 
potential  were  derived  using  a  multi-structure-fitting  procedure,  described  elsewhere 
[18,  19],  and  are  reported  in  Table  1. 
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Table  1.  Potential  parameters  used  [19, 20]. 


Species 

A(eV) 

P(A) 

C  (eVA6) 

O2'  -  O2' 

9547.96 

0.2192 

32 

Al3+  -  O2’ 

1365.79 

0.30096 

2.538 

Cr3+  -  02‘ 

1452.28 

0.30918 

4.472 

Ga3+  -  O2' 

1456.72 

0.30988 

4.616 

Fe3+  -  02‘ 

1478.98 

0.31306 

6.960 

Sc3+  -  O2' 

1587.95 

0.32190 

8.143 

In3+  -  02‘ 

1595.65 

0.32960 

7.402 

Lu3+  -  O2' 

1561.36 

0.33854 

10.01 

Yb3+  -  O2' 

1616.68 

0.33798 

13.34 

Er3+  -  O2' 

1678.21 

0.33781 

10.34 

Y3+  -  O2' 

1721.23 

0.33821 

10.29 

Gd3+  -  O2' 

1868.75 

0.33880 

11.62 

Eu3+  -  02‘ 

1886.71 

0.33975 

11.997 

Sm3+  -  O2- 

1944.44 

0.34080 

12.49 

Nd3+  -  O2' 

1979.11 

0.34148 

13.07 

Pr3+  -  02‘ 

2025.54 

0.34270 

13.83 

Ce3+  -  O2' 

2034.18 

0.34380 

15.86 

Tb3+  -  O2' 

1818.00 

0.33845 

14.33 

Dy3+  -  O2' 

1767.64 

0.33770 

10.94 

Ho3+  -  O2' 

1726.29 

0.33776 

10.72 

La3+  -  O2' 

2051.32 

0.34585 

15.51 

Pu3+  -  O2’ 

1152.00 

0.37480 

12.10 

Table  2.  Shell  Parameter  for  O2  [19, 20]. 


Species 

Y  (e) 

k  (eVA'2) 

O2- 

-2.04 

6.3 
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In  order  to  investigate  the  radiation  tolerance  of  the  perovskite  materials,  lattice 
defects  must  be  considered,  and  their  relative  energies  assessed.  However,  defect 
energies  involve  a  contribution  due  to  lattice  relaxation.  This  relaxation  is  greatest 
close  to  the  defect.  Therefore,  the  lattice  is  partitioned  into  spherical,  concentric 
regions  [21].  The  defect  is  placed  at  the  centre  of  Region  I  and  the  ions  are  relaxed 
explicitly  using  a  Newton-Raphson  procedure.  Region  Ha  is  an  interfacial  region  in 
which  the  forces  between  ions  are  determined  via  the  Mott-Littleton  approximation 
[21,  22]  and  all  ions  are  relaxed  in  a  single  step,  whilst  the  interaction  energies 
between  the  ions  in  region  Ila  and  region  I  are  calculated  explicitly.  Finally,  the  outer 
region  lib  is  effectively  a  point  charge  array  that  provides  the  Madelung  field  of  the 
remaining  crystal.  The  relaxation  energy  of  region  lib  is  approximated  via  the  Mott- 
Littleton  approximation  [23].  In  these  calculations,  the  radius  of  Region  I  is  10  A  and 
that  of  Region  Ila  is  31  A.  All  calculations  were  carried  out  using  the  GULP  code  [24, 
25]. 

The  electronic  polarisability  of  oxygen  ions  is  accounted  for  via  the  shell  model  [26]. 
For  the  purposes  of  this  study  we  have  used  a  massless  shell  with  charge  Y  =  -2.04  |e|, 
is  coupled  to  a  core  of  charge  X=+0.04  |e|,  via  an  isotropic  harmonic  spring  with  force 
constant,  k  =  6.3  eVA'2. 
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Results  and  Discussion 


Crystallography 

Our  traditional  view  of  the  perovskite  lattice  is  that  it  consists  of  small  B  cations 
within  oxygen  octahedra,  and  larger  A  cations  which  are  XII  fold  coordinated  by 
oxygen.  This  structural  family  is  named  after  the  mineral  CaTi03  which  exhibits  an 
orthorhombic  structure  with  space  group  Pnma  [7].  For  the  A3+B3+C>3  perovskites  the 
most  symmetric  structure  observed  is  rhombohedral  R3c  (e.g.  LaA103)  which 
involves  a  rotation  of  the  B06  octahedra  with  respect  to  the  cubic  structure.  However, 
this  distortion  from  the  perfect  cubic  symmetry  is  slight  [27]. 

One  measure  of  the  propensity  for  a  composition  to  adopt  a  high  symmetry  structure 
is  the  classical  tolerance  factor,  t.  This  is  defined  in  equation  2  where  Ra,  Rb,  and  Ro 
are  the  relative  ionic  radii  of  the  A  site  cations,  B  site  cations,  and  oxygen  ions 
respectively  [28]. 


I  _  +Ro 

V2  (Rb  +  Rq) 


Equation  2 


The  tolerance  factor  will  be  equal  to  1  if  the  A  and  B  cations  fit  exactly  into  a  perfect 
cubic  packing  of  oxygen  octahedra.  It  is  therefore  not  surprising  that  LaAlC>3  adopts 
the  highly  symmetric  R3c  structure  since  its  tolerance  factor  is  1.00,  based  on  the 
appropriate  VI  and  XII  coordinate  radii  of  Shannon  [29].  As  the  A  cation  radius 
decreases  and/or  B  cation  radius  increases,  the  tolerance  factor  decreases.  In  the 
perovskite  structure,  this  is  associated  with  the  octahedral  tilting  to  yield  lower 
symmetry  arrangements  which,  in  the  first  instance,  gives  rise  to  an  orthorhombic 
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structure  with  space  group  Pnma.  Even  greater  deviations  lead  to  a  structure  with 
hexagonal  P6scm  crystallography  [30].  In  this  final  so  called  perovskite  variant,  the 
distortions  are  so  great  that  the  A  cations  are  now  essentially  VII  coordinate  and  the  B 
cations  are  V  coordinate.  As  such  it  is  not  clear  that  this  so  called  perovskite  is  very 
perovskite  like! 

Beyond  this  hexagonal  regime,  all  resemblance  to  the  perovskite  structure  is  lost  and  a 
cubic  bixbyite  is  formed  with  Ia3  symmetry.  Here  the  A  and  B  cation  radii  are  very 
similar,  both  cation  sites  are  VI  coordinated  by  oxygen.  Although  dissimilar  to 
perovskite,  it  is  possible  to  consider  bixbyite  in  terms  of  a  defective  fluorite  structure 
[31]- 

In  order  to  elucidate  trends,  we  must  compare  all  compositions  in  a  single  figure.  This 
can  be  achieved  if  the  A  cations  are  ordered  by  size  (assuming  the  XII  coordinate  radii 
values  throughout)  along  an  x-axis  and  the  B  cations  (using  the  VI  coordinate  values) 
along  a  y-axis,  thereby  forming  a  2D  grid  of  compounds.  Such  a  plot  is  shown  in 
figure  1  for  all  96  compounds  for  which  calculations  were  performed.  Of  course,  for 
the  P63cm  and  Ia3  structures,  the  cations  are  not  XII  and  VI  coordinated,  however  to 
allow  for  comparison  the  XII  and  VI  coordinate  radii  are  used  throughout.  In  this  case 
therefore,  the  radii  simply  represent  an  order  parameter. 

The  lattice  energy  for  each  compound  was  calculated  according  to  equation  1, 
assuming  each  of  17  possible  crystal  structures  reported  in  the  literature  for  materials 
with  the  B03  stoichiometry  (including  the  four  mentioned  above).  From  these  data  the 
lowest  lattice  energy  is  identified  for  each  compound  so  that  an  energy  value  is 
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assigned  to  each  composition.  Each  composition  is  also  identified  in  figure  1  by  a 
symbol  which  indicates  its  lowest  energy  structure  (e.g.  a  square  for  bixbyite  etc.). 
Lattice  energy  contours  were  then  generated  using  software  [32],  which  connects 
equal  energy  values  over  the  cation  radii  surface.  Such  maps  facilitate  comparisons 
and  allow  identification  of  compositional  ranges  that  are  of  particular  significance.  In 
this  case  it  is  at  once  clear  that  the  compounds  fall  into  four  groups:  rhombohedral 
R3c,  orthorhombic  Pnma,  hexagonal  P63cm,  and  bixbyite  Ia3.  In  this  regard  it  is 
important  to  note  that  although  17  different  crystal  structures  have  been  reported  for 
ABO3  compounds  in  general,  some  have  only  been  reported  for  A2+B4+C>3  materials 
and  others  represent  only  slight  distortions  from  those  already  mentioned. 

Our  predictions  of  structural  types  are  in  overall  agreement  with  experimental  data 
and  are  consistent  with  changes  in  the  tolerance  factor  (see  figure  2).  In  particular, 
where  we  predict  an  orthorhombic  structure,  experimental  structures  are  known  to  be 
orthorhombic  in  most  cases.  The  same  holds  true  for  the  hexagonal  and  rhombohedral 
structures.  However,  within  the  orthorhombic  and  hexagonal  regimes,  certain 
compounds  are  reported  experimentally  to  exhibit  slight  variations  from  Pnma  and 
P63cm  respectively.  Furthermore,  a  few  compositions  in  the  centre  of  the  map,  such  as 
YFeC>3,  do  not  form  an  ABO3  material,  but  disproportionate  into  garnet  structures. 

Disorder  Processes 

Once  the  lowest  energy  lattice  structure  has  been  identified,  the  individual  component 
defect  energies  were  calculated.  Then  each  compound  can  be  assigned  Schottky, 
oxygen  Frenkel  or  antisite  defect  process  energies  by  following  the  total  defect 
processes  described  in  equations  3,  4  and  5  below. 
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Equation  3 


”nuir->  V"'  +V'"+3V‘‘ 


Ox0  -»  Of  +  Vq  Equation  4 

A*  ->  A‘“  +  Vf  Equation  5 

The  calculated  values  for  A  and  B  cation  Frenkel  energies  are  far  higher  than  those  for 
Schottky,  oxygen  Frenkel  and  antisite  reactions,  and  as  such  will  not  be  considered 
further.  In  a  similar  manner  as  described  for  perfect  lattice  energies,  defect  process 
energies  can  be  represented  over  the  full  compositional  range  on  an  energy  contour 
map  (see  figures  3,  4  and  5  which  use  the  same  energy  contour  interval).  Low  energy 
regions  identify  compounds  for  which  the  defect  process  is  more  favourable.  In  order 
to  allow  for  comparison  between  the  different  defect  reactions,  the  energies  presented 
are  normalised  per  defect,  as  dictated  by  a  mass  action  analysis  [33],  i.e.  the  Schottky 
process  has  a  normalisation  factor  of  5,  whilst  the  antisite  and  oxygen  Frenkel 
processes  have  normalisation  factors  of  2. 

Returning  to  the  Schottky,  oxygen  Frenkel  and  antisite  reactions,  comparing  figures  3, 
4  and  5,  it  is  at  once  apparent  that  the  Schottky  and  oxygen  Frenkel  disorder  processes 
show  less  variation  in  energy  as  a  function  of  cation  radii  than  the  antisite  process. 
Although  the  tolerance  factor  was  originally  derived  to  predict  the  distortion  of  a 
perovskite  lattice  from  an  idealised  cubic  structure,  we  can  also  use  it  to  understand 
the  variation  exhibited  by  the  antisite  energy.  This  is  because  the  tolerance  factor  is 
actually  a  relative  measure  of  the  difference  in  size  of  the  A  and  B  cations  and  by 
inference  how  much  distortion  would  be  necessary  to  accommodate  an  A  cation  at  a  B 
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site  and  vice-versa  (i.e.  as  in  the  antisite  reaction).  In  this  regard,  we  find  it  convenient 
to  quantify  this  for  the  perovskite  compositions  alone,  such  that  low  t  values  (i.e.  less 
than  0.8)  indicate  that  the  cations  fit  fairly  well  in  each  others  sites,  while  high  t 
values  (i.e.  greater  than  0.8)  suggest  that  they  fit  rather  less  well.  Thus  the  value  for  a 
perfect  cubic  perovskite  (i.e.  1)  would  imply  a  high  antisite  energy.  By  comparing 
figures  2  and  5,  it  is  clear  that  there  is  an  excellent  correlation  between  our  calculated 
antisite  energy  and  the  tolerance  factor. 

The  smaller  predicted  variation  in  Schottky  energy  as  a  function  of  cation  radius  is  a 
consequence  of  competing  factors.  In  the  Schottky  reaction  (equation  3)  the  displaced 
ions  form  new  lattice  and  as  such,  the  lattice  energy  and  hence  its  variation  with 
component  cation  radii  is  an  important  factor.  This  was  shown  in  figure  1  and  clearly 
the  lattice  energy  gain  is  greatest  for  compositions  in  the  bottom  left  hand  comer.  The 
second  term  concerns  the  energy  to  form  vacancies.  Higher  energies  (less  favourable) 
will  result  in  materials  with  higher  Madelung  or  electrostatic  energies.  Although  these 
vary  somewhat  between  structure  types,  values  principally  increase  as  the  interatomic 
separation  decreases.  As  such  we  expect  higher  energies  in  the  bottom  left  hand 
comer,  i.e.  the  opposite  variation  observed  with  the  lattice  energy  (see  figure  1). 
These  electrostatic  terms  are  therefore  working  in  opposition,  and  the  net  effect  is  a 
smaller  variation  in  Schottky  energy  as  a  function  of  cation  radii.  Contributions  from 
the  short  range  pair  energies  now  become  important.  Consequently  structure  type 
boundaries  where  there  is  a  change  in  ion  coordination  should  show  up  in  the 
Schottky  energy  map.  In  figure  3,  we  show  evidence  of  such  a  change  between  Pnma 
(VI,  XII)  and  P63cm  (V,  VII)  but  not  R3 c  (VI,  XII)  to  Pnma (VI,  XII). 


11 


By  comparing  the  results  in  figures  3  and  5  it  is  at  once  evident  that  there  is  a  change 
in  lowest  disorder  process  energy  from  Schottky  to  antisite  as  a  function  of  the  cation 
radii.  If  we  express  this  in  terms  of  the  tolerance  factor,  the  swap-over  occurs  at  a  t 
value  between  0.99  and  0.98  with  the  antisite  reaction  dominant  for  all  lower  values. 
Thus,  only  in  a  small  portion  of  the  cation  radii  space  (i.e.  at  the  bottom  right)  does 
Schottky  disorder  dominate  over  antisite.  Furthermore,  in  this  region,  the  difference  in 
process  energies  is  not  great,  and  thus  materials  such  as  LaAlCb  should  exhibit 
complex  defect  equilibria.  We  can  also  compare  the  oxygen  Frenkel  and  antisite 
energies.  The  oxygen  Frenkel  energy  decreases  somewhat  as  the  A  cation  radius 
increases  and  the  B  cation  radius  decreases,  in  a  similar  manner  to  the  antisite  process 
energies.  The  oxygen  Frenkel  energy  decreases  somewhat  as  the  A  cation  radius 
decreases  and  the  B  cation  radius  increases,  in  a  manner  similar  to  the  antisite  process. 
Although  it  does  drop  to  a  value  considerably  below  the  Schottky  for  a  majority  of  the 
compositional  range,  the  oxygen  Frenkel  energy  does  not  decrease  to  the  same  extent 
as  the  antisite  and  is  never  the  dominant  defect  process.  Again  we  see  some  evidence 
of  the  Pnma  to  P63CIT1  structure  change. 

Comparison  to  Pyrochlore/Fluorite 

The  perovskite  defect  process  energies  can  now  be  compared  to  the  equivalent  lowest 
energy  defect  process  energies  in  other  materials.  Here  we  can  make  that  comparison 
to  the  fluorite  and  pyrochlore  oxides  investigated  previously  [12,  19].  In  this  case,  the 
antisite  energies  of  the  perovskites  are  well  above  those  for  the  zirconium  containing 
pyrochlore/fluorite  materials  (e.g.  E^Z^O?  (0.6  eV)  and  Gd2Zr2C>7  (1.8  eV)  [12]).  As 
such,  only  selected  P6acm  perovskites  will  exhibit  defect  concentrations  as  high  as  the 
zirconate  pyrochlores.  On  the  other  hand,  the  titanate  pyrochlores  exhibit  antisite 
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defect  energies  of  typically  2.9  eV  which  is  equivalent  to  the  lowest  energy  for  the 
Pnma  perovskites. 

The  oxygen  Frenkel  defect  energies  for  the  perovskites  (figure  4)  are  again  higher 
than  those  for  the  pyrochlores.  For  example,  the  zirconate  pyrochlores  Er2Zr207, 
Gd2Zr2C>7  exhibit  oxygen  Frenkel  energies  of  2.5  eV  and  2.0  eV  respectively  [19], 
both  being  lower  than  any  of  the  perovskite  materials.  However,  the  titanate 
pyrochlores  have  values  that  are  similar  to  the  Pnma  perovskites,  for  example, 
Gd2Ti207  has  an  oxygen  Frenkel  energy  of  5.6  eV  [19].  Consequently,  values  for  the 
bixbyites,  which  are  as  low  as  3  eV,  compare  favorably  to  the  titanate  pyrochlores  and 
are  not  significantly  higher  than  the  zirconates. 

Up  to  this  point  we  have  suggested  a  model  that  considered  only  isolated  defects. 
However,  the  stability  of  such  defects  depend  not  only  on  their  individual  energies, 
but  also  on  the  interactions  between  them.  This  is  characterised  by  the  reduction  in 
formation  energy  when  the  defects  are  associated  (known  as  the  binding  energy).  The 
work  on  pyrochlores  considered  the  representative  defect  cluster 
{Ag  :B*  \  V"  :0"}x ,  a  combination  of  antisite  and  oxygen  Frenkel  reactions  where 

the  defects  are  in  their  nearest  neighbour  positions  [19].  Interestingly,  when  different 
geometries  of  this  cluster  are  investigated  for  pyrochlore,  the  oxygen  Frenkel 
component  somtimes  self  annihilates  around  the  cation  antisite  pair.  Such  behaviour 
was  not  observed  in  equivalent  perovskite  clusters.  This  may  in  itself  be  a  significant 
difference  between  pyrochlores  and  perovskites  in  regard  to  their  ability  to 
accommodate  large  defect  populations.  Nevertheless,  in  order  to  compare  equivalent 
defect  clusters,  here  we  simply  compare  energies  for  adjacent  antisite  pairs 
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{Ag  :BX}X  (see  table  3  for  examples).  Critically,  for  perovskites  the  value  when  the 
defects  are  associated  (i.e.  clustered)  are  only  modestly  less  than  when  the  defects  are 
isolated). 

Table  3.  Comparison  of  the  total  defect  energies  for  isolated  defects  and  clusters 
for  different  compositions  of  perovskite  and  pyrochlore  oxides.  N.B.  the  energies 
in  this  table  have  not  been  normalised  to  facilitate  comparisons. 


Total  Defect  Energy  (eV) 

{Axb 

:  B*)X 

Isolated 

Clustered 

Perovskite 

GdAlQ, 

8.20 

7.05 

LaAlOs 

11.98 

10.59 

GdFe03 

4.24 

4.03 

LaFeO;* 

7.44 

6.30 

Pyrochlore 

[2] 

Gd2Ti207 

6.09 

4.98 

LftaTiaOr 

5.94 

4.59 

Gd3Zr30T 

3.68 

2.86 

La2Zr207 

4.52 

3.45 

Conversely,  in  pyrochlore  materials,  cluster  formation  values  are  considerably  less 
than  for  isolated  defects  (again  see  table  3).  The  larger  binding  energy  for  pyrochlores 
now  means  that  we  see  a  significant  difference  between  the  antisite  energy  for 
GdFeCh  compared  to  lower  zirconate  pyrochlore  values.  The  other  perovskites 
maintain  their  higher  defect  process  energies. 


Clearly,  for  the  pyrochlore  materials,  these  antisite  defects  are  more  strongly  bound 
than  for  the  perovskite  materials.  This  may  well  imply  that  in  perovskites,  as  the 
defect  concentrations  increase,  the  retained  energy  increases  more  quickly.  On  this 
basis,  all  the  perovskites  should  accommodate  lower  defect  concentrations  than 
pyrochlores.  However,  to  be  definitive  would  require  better  analysis  probably  using  a 
molecular  dynamics  approach  (because  of  the  great  number  of  cluster  configurations). 
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Conclusion 


As  a  consequence  of  the  broad  compositional  range  of  ABO3  materials  considered,  the 
assumed  crystal  structures  of  individual  compositions  varies  considerably  from 
rhombohedral  R  3  c  when  the  A  and  B  cation  radii  are  most  different,  to  orthorhombic 
Pnma,  to  hexagonal  P6scm,  and  finally  to  the  cubic  bixbyite  Ia3  as  the  constituent 
cation  radii  converge.  Despite  this  drastic  change  in  the  crystallography  of  the 
system,  a  smooth  transition  in  the  lattice  energy  across  these  structural  boundaries  is 
predicted  (see  figure  1). 

As  a  first  approximation  to  understanding  the  variation  in  crystal  structure,  the 
classical  tolerance  factor  provides  an  important  benchmark.  The  tolerance  factor 
follows  the  changes  remarkably  well  and  infers  that  the  atomistic  simulation  results 
are  in  agreement  with  classical  perovskite  theory  (see  figure  2). 

The  structural  changes  also  manifest  themselves  in  the  defect  chemistry  of  these 
ABO3  materials.  The  defect  process  energies  of  the  Schottky  and  oxygen  Frenkel 
reactions  show  abrupt  changes  in  compositional  regions  (i.e.  cation  radii  space) 
associated  with  changes  in  crystal  structure  (see  figures  3  and  4).  Conversely  in  the 
antisite  reaction  energy  decreases  smoothly  and  rapidly  as  a  function  of  the 
component  cation  radii  across  crystal  structure  boundaries  (see  figure  5). 

For  compounds  with  large  A  cation  radii  and  small  B  cation  radii  (i.e.  the  lower  right 
hand  comer  of  the  maps)  the  Schottky  defect  reaction  is  the  dominant  (lowest  energy) 
process,  whereas  for  the  rest  of  the  compositional  range  the  antisite  process 
dominates.  However,  if  one  considers  defect  processes  on  just  the  oxygen  sublattice, 
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there  is  also  a  change  over  in  the  dominant  reaction  from  Schottky  to  oxygen  Frenkel 
that  occurs  close  to  the  orthorhombic  Pnma  to  hexagonal  P63cm  crystal  structure 
boundary.  Such  process  mechanisms  are  important  because  such  defect  reactions  on 
the  oxygen  sublattice  control  the  equlibria  which  provide  the  defects  that  mediate 
transport  mechanisms  of  oxygen  through  the  lattice. 

In  an  attempt  to  understand  the  interplay  between  defects  on  the  cation  and  anion 
sublattices,  defect  clusters  have  been  considered  that  combine  a  cation  antisite  defect 
pair  with  an  oxygen  Frenkel  pair.  The  clusters  in  perovskite  have  been  compared  to 
equivalent  clusters  investigated  previously  in  the  pyrochlore  materials.  One  of  the 
startling  differences  between  the  two  structures  is  that  within  the  pyrochlore  materials, 
different  geometries  of  the  cluster  lead  to  the  annihilation  of  the  oxygen  Frenkel  pair 
vastly  reducing  the  energy  of  the  cluster.  Such  mechanisms  were  not  seen  in  the 
perovskite  materials.  Partly  as  a  consequence,  a  different  cluster  formed  from  an 
adjacent  antisite  defect  pair  was  considered  to  allow  for  comparison  between  the 
perovskite  and  pyrochlore.  Here,  the  energies  for  clustered  defects  in  the  perovskites 
are  only  modestly  lower  than  when  the  defects  are  isolated,  conversely,  the  equivalent 
cluster  within  the  pyrochlore  materials  is  much  lower  in  energy  than  the  isolated 
defects.  As  such  the  antisite  defects  are  more  strongly  bound  together  in  the 
pyrochlore  material  than  the  perovskite.  This  implies  that  as  the  defect  concentrations 
increase  within  perovskites,  the  retained  energy  increases  rapidly  and  as  such  the 
perovskites  should  not  be  able  to  accommodate  such  large  defect  concentrations  as  the 
pyrochlore  material 
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Figure  2.  Classical  Tolerance  Factor  variations  asa  function  of  cation  radii. 
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Figure  3.  Schottky  reaction  energy  Contour  map. 
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Figure  4.  Oxygen  Frenkel  reaction  energy  contour  map. 
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Figure  5.  Antisite  reaction  energy  contour  map. 
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